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We use lattice techniques to study the closed-string spectrum of SU (N) gauge theories in 2 + 1 
dimensions. We calculate the energies of the lowest lying ~ 30 states for strings with lengths 
Q-ij between I ~ 0.45 fm and I ~ 3 fm, and compare to different theoretical predictions. We obtain 

, ' unambiguous evidence that the closed-strings are in the universality class of the Nambu-Goto free 

bosonic string. Moreover, we clearly see that our data can be described by a covariant string theory 

>. 

with a small/moderate correction down to very short distance scales, and possibly on all distance 
scales at large- N. 
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I. INTRODUCTION 



The generation of electric flux tubes in confining gauge theories is a basic phenomena 
that characterises the vacuum of these theories. In this paper we study the energy spectrum 
of these flux tubes using lattice techniques. We are mainly motivated by the following two 
reasons. 

Firstly, a flux-tube whose length I is much larger than its width, is expected to be 
described by an effective low energy string action S fj. Establishing what is the structure 
of S'efr is of fundamental interest, and can be done by studying the energy spectrum of the 
flux-tube. More precisely, while the energy levels E n of a string with tension a obey 

E n (l) l =°°al, (1.1) 

for finite I there are corrections to Eq. (11. ip that reflect various properties of S'efr- I n par- 
ticular, the universal coefficient of the 0(1/1) 'Luscher-term' is determined by the number 
of the massless degrees of freedom propagating on the worldsheet of the string, thus re- 
vealing the IR universality class of S e g In contrast, other details of the spectrum, like 
its degeneracies and the form of the subleading contributions to the Luscher term, are in 
general not universal and depend on the particulars of the non-renormalizable terms in S'efr- 
Consequently, by studying the /-dependence of the string spectrum, one can directly learn 
about S e ff, and about the length String above which S e fr begins to be a good description and 
a string is formed. 

Our interest in the flux-tube spectrum is also driven by the following more practical 
reason. An essential step in any lattice study is to calculate the lattice spacing a in physical 
units. This is sometimes done by extracting the dimensionless combination a^fo from a 
flux-tube's ground state energy given in lattice units aE n=Q . Consequently, this requires 
that we know how Eq depends on a and I. To date, this dependence was approximated by 
correcting Eq. (11.11) with just the Luscher term. For the purpose of high precision lattice 
studies, however, this approximation can become insufficient and it is imperative to know 
what are the subleading corrections to Eq. (11.11) that go beyond the Luscher term. 

In this work we focus on the spectrum of closed flux-tubes in pure SU(N) gauge theories 
in D = 2 + 1 dimensions. The flux-tubes that we study have lengths that range from 
/ ~ 0.45fm to I ~ 3fm, while the lattices we use have spacings that range from ~ 0.06fm to 
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~ 0.2fm, depending on the value of N and I. (Here, despite working in 2 + 1 dimensions and 
in pure gauge theories, we choose to define lfm through the convention a = (440 MeV) 2 .) 

The gauge groups that we study have N = 2,3, 4, 5, 6, 8, and here we are motivated by 
the special role that the large- iV limit plays in the physics of confinement; From a field 
theory point of view, several processes that cannot be described by a simple low energy 
effective string theory, such as glueball-string mixing, string/anti-string mixing and de- 
confinement /instability of short strings, become less important with increasing N (see the 
discussion in Section [II C 3|) . This is expected to make the description of the flux-tube in 
terms of S e s better and simplify the form of the latter. From a more general string/gauge 
duality point of view (for example see 

SB) 

, it is the large- N limit of QCD that one may 
hope to describe by a string theory. This makes non-perturbative information on this limit 
important both to guide the searches for such a dual string theory, as well as to understand 
it beyond the supergravity approximation. 

There are two stages to our calculation which henceforth we refer to as A and B. In stage 
A we perform high precision measurements of the energy of the flux-tube's ground state. In 
particular, we aim to isolate the different corrections to Eq. ( II. ip and, while controlling them, 
extract a value for a. We then use the latter to compare our data to different theoretical 
predictions, and especially to the Nambu-Goto (NG) free string model. This stage of our 
work has two immediate practical implications. Firstly, it allows us to test the analytic 
work in [4| which predicts that in D = 2 + 1 the 1/7 3 correction to E n has a universal 
coefficient, being precisely that predicted by the NG model. 1 Secondly, a precise control of 
the corrections to Eq. (11. ip has enabled us to perform, in a companion study j^K a precise 
test of the Karabali-Kim-Nair prediction for the value of a in 2 + 1 dimensions 6(. 

In stage B we put our emphasis on the excites states in the spectrum and measure the 
energy of the lowest ~ 30 states. This allows to extend the comparison with theoretical 
predictions to states with more quantum numbers and a nontrivial degeneracy structure. 

The study of confining flux-tubes with lattice techniques has been an active field of 
research for the past three decades, and we refer the reader to some recent papers Q. 
These include works that vary in the selection of the gauge group, the number of space-time 



1 Note that, using the method of |8|], the authors of [9( assert that this universality persists for a general 
value of D. 
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dimensions, and the boundary conditions imposed on the flux-tube. Other papers which are 
particularly relevant to our current study are mentioned later on in the text. 

The following is the outline of this letter. We begin in Section [III by describing our 
methodology, proceed to Section [TTT] where we discuss the theoretical expectations for the 
string spectrum, and move to Sections IIV11 VII where we present the results of the calculations. 
In Section IVIII we summarise our results and make a few remarks on their theoretical and 
practical implications. 



II. METHODOLOGY 



In this section we describe a few aspects of our methodology. We begin with the lattice 
construction, proceed to discuss our general strategy to extract the energy spectrum from 
correlation functions, and end by listing the main systematic errors and by expanding on 
how we control them. 



A. Lattice construction 



We define the gauge theory on a discretized periodic Euclidean three dimensional space- 
time lattice, with spacing a. The fields are SU(N) matrices assigned to the links of the 
lattice, and the Euclidean path integral is given by 

Z = J DUexp{-(3S w ). (2.1) 

Here j3 is the dimensionless lattice coupling, and for our action is related to the dimensionful 
coupling g 2 by 



2N 

lim/3=— . (2.2) 
a->(r ag 2 v ' 

In the large- N limit, the 't Hooft coupling A = g 2 N is kept fixed, and so we must scale 

(3 = 2N 2 /X oc iV 2 in order to keep the lattice spacing fixed (up to 0(1/N 2 ) corrections). 

The action we choose to use is the standard Wilson action 



p 



^-ReTr U P 

N 



(2.3) 



where P is a lattice plaquette index, and Up is the plaquette variable obtained by multiplying 
link variables along the circumference of a fundamental plaquette. We calculate observables 



by performing Monte-Carlo simulations of Eq. fl2.ll) . in which we use Cabibbo-Marinari 
updates of the link matrices with a mixture of Kennedy-Pendelton heat bath and over- 
relaxation steps for all the SU(2) subgroups of SU(N). 

In stage A we have studied N = 2,3,4,5,6,8 with three lattice spacings, a ~ 
0.06,0.08,0.11 fm. In stage B we studied SU(3) with a ~ 0.04,0.08 fm, and SU{6) with 
a ~ 0.08 fm. The string lengths / in both stages ranged between ~ 0.45 fm and ~ 3 fm, 
depending on the values of N and a. For more details on the lattice parameters of our field 
configurations see Tables [THTTl 



B. General strategy 

Since we are mainly interested in the way the flux-tube spectrum reflects the properties 
of a standard low energy string theory of the NG type, we restrict ourselves to the spectrum 
of closed flux-tubes. By this restriction we avoid a class of short- distance contributions 
to the energies, such as the Coulomb interaction between sources, that cannot be easily 
accommodated in an effective string theory (see the discussion in Section III C 3p . 



We calculate the energies of flux tubes that are closed around a spatial torus. We do so 
from the correlators of suitably smeared Polyakov loops that wind around one of the spatial 



tori and that have vanishing transverse momentum. This is a standard technique [10J, 111 ] 
with the smearing/blocking designed to enhance the projection of our operators onto the 
ohysical states. (We use a scheme that is the obvious dimensional reduction of the one in 



12j.) We classify our operators using the following quantum numbers 

1. P — ± : Parity in the direction transverse to Euclidean time and to the string contour. 

2. q = 0, ±1, ±2, . . . : The momentum, in units of 2ir/l, along the string. 

For each combination of these quanta we construct the full correlation matrix over our space 
of loop operators, and use it to obtain best estimates for the string states using a variational 

* - nnnn 

method applied to the transfer matrix T = e aH - again a standard technique 10, llll, ll2|, ll4| . 

The size iV of our N Q x iV correlation matrices depend on the states we are interested 
in. In stage A, where we focus on the lowest state, our correlations are built from Polyakov 
loops that wind once around the torus in straight paths and in all possible blocking levels. 
Consequently, this means that we restrict ourselves to states with q = and P = +, and this 



provides us with N Q = 3 — 5 states. In stage B we include Polyakov loops with many different 
paths (wave-like paths, pulse-like paths, etc.). This increases our number of operators to 
N D ~ 80 — 200 and allows us to probe states with negative P and nonzero q. Finally, an 
extension of our calculation for the ground state to w > 1 will be reported in [15I ]. 

Once we obtain the string energies for different values of /, we fit the ground state energy 
in powers of l/l and obtain an estimate for the function E (l, a). With this empirical formula 
we extract the string tension in lattice units, a 2 cx, for each data set. Substituting the results 
in different theoretical predictions for the spectrum, we conclude by comparing the latter 
with our data. 



C. Systematic errors 

We control several systematic errors in different stages of our study. These include the 
way we extract the string energies from correlation functions, the way our results approach 
the infinite volume and continuum limits, and the way one may interpret the results for 
small values of I as reflecting an effective string picture. In the next three sub-sections we 
expand on each systematic error, and on the way that we control them. 



1. Extracting string energies from correlation functions 

The output of the variational technique is a set of operators that couples "best" to a 
set of low energy states. For example, our best operator for the string's ground state has 
typically an overlap ~ 99% onto that state so that the normalised 'ground state' correlation 
function satisfies 

C(t) = (l-\e\)exp{-E (l)t} + \e 1 \exp{-E 1 (l)t} + ... ; £ = \e\ ~ 0.01 (2.4) 

where E , E\ are the ground and first excited state string energies. (Since our time-torus is 
finite, we use cosh fits rather than simple exponentials, although in practice we use L t large 
enough for any contributions around the 'back' of the torus to be negligible.) To extract E 
from this correlator one can fit with a single exponential for t > t , discarding t < to, and 
choosing to to be the minimum value so that a statistically acceptable fit is obtained. This is 



a reasonable approach and one followed in lOj, QJ|. However it neglects the systematic error 
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arising from the fact that there is certainly some excited state contribution as demonstrated, 
for example, by the fact that one cannot obtain a good fit with a single exponential from 
t = 0. To control this systematic error we also perform fits with two exponentials, with 
a fixed mass E* for the excited state, resulting in a mass Eq(E*) for the ground state. 
Typically Eq{E*) is smallest when E* is as small as possible, i.e. E* = E%, and is largest 
when E* = oo, i.e. effectively a single-exponential fit. So the true value typically satisfies: 



From here on, we refer to the single-cosh fitting procedure by 'S', and to the double-cosh 
fitting procedure, by 'D'. This particular systematic error becomes more important as the 
string length I increases because the overlap of our lattice operators typically decreases with 
increasing I, and the excited state energy approaches the ground state energy. In practice, 
we find that as long as I < 5/ ' \fa one can neglect this systematic error at the level of our 
current statistical errors. Consequently we present results from the 'D' fitting procedure 
only when I > 5/ \fa. 

2. Finite volume and discretisation effects 

To avoid finite volume effects in the calculation of the closed string spectrum we follow 
jl^ and increase the transverse and temporal extents of the torus, when we decreases the 
length of the string /. We performed explicit finite volume checks for a restricted set of 
parameters, and found that the transverse volumes used in our study are large enough to 
avoid any observable effects at the level of our statistical accuracy. To check for finite lattice 
spacing effects we perform several of our calculations with different lattice spacings. 

3. A string interpretation of the flux-tube spectrum for short lengths ? 

The energy spectrum of confining flux-tube is expected to be much more complicated than 
that of a simple effective string and to approach the latter only at large /. Therefore, in an 
ideal calculation one would study the l/l terms in the string energies only for / ^> \j\fa. 
In practice, however, the l/l terms are numerically small and a useful measurement at large 
/ would require an unrealistically large statistical sample. Instead, we study strings with 
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l/x/o" < / < 6/a/o" for which the corrections are not negligible and can be reasonably fitted. 
This procedure is not free of ambiguities since the short-Z spectrum may be sensitive to 
phenomena not accommodated in an effective string theory, but as we explain below, it is 
unlikely that this ambiguity is significant in our calculation. 

The first obvious source of short-distance "contamination" in studies of flux-tube spectra 
is the presence of a Coulomb interaction between the static charges at the ends of open flux- 
tubes. This can be a significant portion of the total energy when the flux tube is short, so it 
is not clear whether the deviations from the infinite length limit seen in the open channel (for 
example see [l7] and reference within) are due to this Coulomb interaction or whether they 
reflect string interactions. Our calculation, however, is free of this contamination simply 
because we study closed flux-tubes and by construction these do not have static charges and 
a Coulomb interaction. 

A different type of non-stringy phenomenon that does occur in the closed channel is the 
deconfinement transition. In 2 + 1 dimensions this finite temperature transition takes place 

n 

at a temperature T = Td — 0.9/ y/a [18(. By identifying the length of the compact direction, 
I, with the inverse temperature T _1 , this means that our calculation necessarily breaks down 
when I < 1/Td — 1.11/a/cx, but also that interpreting the flux-tube energies as coming from 
the dynamics of S e g may be questionable in the confined phase, when I ~ (1/T^) + . This 
is particularly true for N = 2, 3 when the transition is second order and the I — > (1/T^) + 
behaviour of Eq(1) will be governed by appropriate critical exponents, and may also occur 
for N = 4,5, where the first order transition is relatively weak. In contrast, for N > 6, the 
transition is strongly first order and it is quite possible that a confining string description 
exists even when / ~ (1/T<i) + . 2 

Another short-distance phenomenon involves glueballs : by decreasing I the energy needed 
to excite the string grows (see Section IHip and if we take the spectrum of the free bosonic 
string prediction as a guide, then the threshold for the first excited state to emit the lightest 
++ glueball and decay to the ground state is reached when I ~ 1.53/ y/a. In practice, 
however, since the amplitudes of these mixing processes are subleading in 1/N 2 , they may 
be suppressed even for SU(3). 



2 As an extreme example, note that for very large values of N and in 3 + 1 dimensions, it is possible to 
study confining flux-tubes even below I/T4 [i^ . 
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To conclude, the short-distance non-stringy phenomena that may contaminate the string 
picture interpretation of the closed flux-tube spectrum, go away at large-iV. Hence, by 
performing our calculations for increasingly large values of N we explicitly check how large 
are these effects and thereby control them. 



III. THEORETICAL EXPECTATIONS FOR E n (L) 

Let us now discuss the theoretical predictions to which we compare the measured flux- 
tube spectrum. 



A. The spectrum of the Nambu-Goto model 



The action of the NG model [20| is the area of the worldsheet swept by the propagation 

is quantum-mechanically consistent only 



211]). but since this anomaly is suppressed 



of the string. Due to the Weyl anomaly this mode 
in the critical dimension D = 26 (see for example 

rn. 

for long strings [22j it can still be considered as an effective low energy model. 

The single string states can be characterized by the number of times w that the string 
winds around the torus. The spectrum in each case corresponds to the transverse oscillations 
of the worldsheet that correspond to movers that travel clockwise and anti-clockwise along 
the string. Thus the string states are characterised by w, by the occupation number riL(R)(k) 
of left (right) movers that carry energy k, and also by the centre of mass momentum p c . m .. By 
projecting to zero transverse momentum {p c . m .)± we are left only with the momentum along 
the string axis which is quantized in units of 2nq/l with q = 0, ±1, ±2, . . . for a string of 
length I. These quanta are not independent of Ul^r and obey the level matching constraint 3 

N L -N R = qw, (3.1) 

where Num enumerates the momentum contribution of the left (right) movers in a certain 
state as follows 

Nl = Y, E Mk)k, iV*=E E n R {k')k'. (3.2) 

fc>0 n L (fc)>0 fc'>0 n R {k')>0 



This condition constraints the physical states to be invariant under the gauged diff-invariance of the action 



21] , and is effectively momentum conservation. 
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It is customary to characterise the string states as irreducible representations of the 
SO(D — 2) symmetry that rotates the spatial directions transverse to the string axis. In 
our D = 2 + 1 dimensional case, this group becomes the transverse parity P and acts by 
assigning a minus sign for each mover on the worldsheet. As a result the string states are 
eigenvectors of P with eigenvalues 

P = (_i)E^i™i(^)+E J= i n fl( fc P _ (3.3) 

Finally, the energy of a closed-string state with the above quanta is (here we write it for 
a general number of spacetime dimensions D) 

ftp \* f 7 V , a (Nl + Nr D-2 \ , f^g\ 2 

B. Effective string theories 

Since in 2 + 1 dimensions the NG string is at best an effective low-energy string theory, 
it makes sense to generalise it and write the most general form of an effective string action 
S c g consistent with the symmetries of the flux-tube system. This was done some time ago 
for the w = 1 and q = in [1] and the spectrum obtained for a general number of space-time 
dimensions D was 

4.71- / n _ o\ 

E n = al + — (n- — — + 0(1// 2 ), (3.5) 



I V 24 

with n = 0, 1, 2, . . .. Here the second term on the right hand side of Eq. (13.51) is known 
as the Luscher term and is expected to be universal and independent of the particulars of 
IR-irrelevant interactions of the low energy effective string theories. Indeed, it can be easily 
verified that the NG model obeys this universality by expanding the square-root of Eq. ( 13. 4ft 
to leading order in 111. 

n n 

The work [l| was more recently extended in |4|, where the authors established and used 
a certain open-closed string duality of the effective string theory. Using this duality they 
showed that for any number of spacetime dimensions the 0(1 /I 2 ) is absent from Eq. (13.51) . 
and that in D = 2 + 1 the 0(l/l 3 ) has a universal coefficient. Consequently, in 2 + 1 
dimensions, Eq. ( 13. 5p is extended to 

^■=" + t("-s)-|f("-s) >+0 ^- (3 ' 6) 
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In the context of the covariant string description such as that of Section IIII Al Eq. 
implies that 

E n = J(al) 2 + 87ra (n - ^) + 0(1/F), (3.7) 

which is a particularly convenient form since the two first terms under the square root are 
the prediction of the NG model for the energy squared, that we find to be a very good 
approximation (see below). A different approach that also leads to similar conclusions is 
the Polchinski-Strominger effective string theory [8:]. While it originally yielded Eq. (13. 5p . it 
was used recently in js] to extend the analysis to higher powers of l/l and leads to the same 
conclusions as 4J, but for all values of D. 

Motivated by these recent developments, we take our main fitting ansatz for the spectrum 
to be 

Ete = ^El G -aj£^ ; p > 3, (3.8) 

where E^q is the NG prediction given by Eq. (13.41) and where C p are dimensionless coeffi- 
cients that in general can depend on the quantum numbers of the state. 

Let us pause and make the following comment about the relation between Eq. ( 13. 6j) 
and Eq. (13. 7\\ . Consider performing a large-/ expansion of the square-root in Eq. ( 13.71) . 
The result would include not only Eq. (13.61) . but also many additional terms that begin at 
0(1/ 1 5 ), and that come from the expansion of the second term under the square-root. To see 
whether these extra terms are important, we can use our data and compare its deviations 
from Eq. ( 13.61) to its deviations from Eq. ( 13. 7p . If we find that the latter are systematically 
smaller, this would then reflect the naturalness of adding a correction term in a covariant 
way (correcting E\ as in Eq. (13.71) and Eq. (13.81) ) rather than correcting the energy E n itself 
(as Eq. (13. 6p does). From a field-theoretical point of view, such as the one in [4|, it seems 
that it would be hard to get a prediction of the form of Eq. (13 .7p since it can be viewed as 
an implicit resummation of an infinite series of powers of the Luscher term. 



IV. A PRECISE MEASUREMENT OF THE GROUND STATE ENERGY E (a,l) 

We performed high precision calculations of the dependence of the ground state en- 
ergy on the length of the string. The calculations were done for the numbers of colors 
N = 2,3,4,5,6,8. The lengths of the strings were restricted to obey I > with the 
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deconfinement transition temperature given in [18(]. As examples, we present the results 
for iV = 3 and N = 6 with (3 = 14.7172 and (3 = 90.00 respectively, in Fig. Q] 




FIG. 1: The ground state energy per unit length, in lattice units, vs. the string length for SU(3) 
and (3 = 14.7172 (left panel) and for SU(6) and = 90.00 (right panel). All the data points are 
the result of the single exponential fits, except when I > 5/y/a, where we use the results from 
double exponential fits (see Section III C lj) . The blue lines are the results of using the fitting ansatz 
Eq. ()3.8p . Using the values we obtain for the string tensions we plot the NG predictions of Eq. (13.41) 
(magenta lines) and of the Luscher formula in Eq. (|3.5p (black lines). 

Fitting the data with Eq. (13.81) and p = 3 we obtain acceptable fits only when we use the 
'D' data for strings with I > 5/y/a. This indicates that for long strings, for which the first 
excited energy is relatively close to the ground state energy, the control over contamination 
of the excited states is crucial. The results for the other gauge groups and lattice spacings 
are similar to the SU(6) results, and we present the parameters and goodness of fits of the 
form in Eq. (13.81) with p = 3 in Table [fl 

From the table one can also see that the iV-dependence and the lattice spacing dependence 
of the correction C3 is moderate/small. We have also performed fits with p = 4, and obtained 
comparable x 2 > which means that with the current data for the ground state, we cannot 
unambiguously determine the power of the correction term in Eq. ( 13.81) . 

As Fig. [T] demonstrates, the Luscher term formula and the NG formula are indistinguish- 
able from each other and agree with the fit at large values of /. The situation at small 
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Configuration details 


a 2 a 


c 3 


Confidence level 


SU{2), 13 = 5.6 


0.074636(25) 


0.0100(41) 


53% 


SU(S), (3 = 14.7172 


0.068127(47) 


0.1633(112) 


87% 


SU(4), f3 = 28.00 


0.063346(73) 


0.2943(282) 


99% 


SU(4), f3 = 50.00 


0.017184(14) 


0.1886(338) 


10% 


SU(5), f3 = 80.00 


0.016874(12) 


0.0554(139) 


68% 


SU(6), (3 = 59.40 


0.077460(81) 


0.1098(169) 


16% 


SU(6), 13 = 90.00 


0.029601(23) 


0.1163(159) 


88% 


SU{8), \3 = 108.00 


0.075351(141) 


0.1865(575) 


68% 


SU{8), \3 = 192.00 


0.020200(24) 


0.0864(109) 





TABLE I: The lattice parameters (3 and a 2 a and C3 in the fit Eq. (|3.8p . The fits were performed 
for the S data sets except for strings of length I > 5 where we used the D type of data (Single 
and Double exponential fitting of the correlation functions - see Section III C ip . For SU{8) and 
(3 = 192.00 our confidence level is small, and comes from a large scatter of the data points around 
the fit. This indicates that our errors are probably underestimated in this case. 

values of I, however, is completely different. There, the Luscher term is clearly insufficient 
to describe the data, indicating the importance of the subleading corrections to Eq. (13. 5ft for 
I < 3/y/a ~ 1.35 fm. In contrast, the NG prediction is very good, even at the lowest value 
of I. This remarkable fact also explains why we choose Eq. (13. 8ft to be our fitting ansatz, 
rather than using a fit which assumes that the energy itself is a power series in 1/7. 
Finally, we perform the following two type of fits 



Fitl. 
Fit 2. 



E (l,a) 
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alf 



ira 



(2) 
eff 5 



(4.1) 
(4.2) 



to pairs of adjacent values of / in plots of the type of Fig. [TJ 

The effective central charge C e g(l) in both type of fits is expected to approach a value of 
unity for large values of I, but to deviate from 1 when I decreases because of the higher order 
0(1 /I 3 ) terms. The fit in Eq. (14.11) is the one usually performed (see for example 
is the analog of the effective central charge fits performed in the open channel [4 
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2j}and 



24J. In 



contrast, the ansatz in Eq. (14 .2p assumes that most of s deviation from 1 comes from 

13 



the l/l terms that are accommodated in the NG prediction. Indeed, this assumption is 
confirmed by the data and approaches 1 much faster than C^jj , and we present both, 
for the 'S' data sets of SU(3) with (3 = 14.7172 and of SU{6) with (3 = 90.00, in Fig. El 



1 1.5 



C „ from Fit 1 

eff 

C „ from Fit 2 

eff 



2 2.5 3 3.5 4 4.5 



1.5 



o 



0.5 



1 1.5 



C _ from Fit 1 

eff 

C _ from Fit 2 

eff 



2 2.5 3 3.5 4 4.5 

IVa 



(1 2) 

FIG. 2: The effective central charges vs. the string length for the 'S' data points, for SU(3) 

and (3 = 14.7172 in the left panel, and for SU(5) and (3 = 80.00 in the right panel. 

Finally, in our largest statistical sample, which was obtained for SU(2), we find 



C$ = 0.9991(55) at / ~ 0.68 fm. 



(4.3) 



V. TEST OF THE NG PREDICTION FOR THE 1ST EXCITED STATE 

Equipped with values for the string tensions, we are now in a position to test the 
parameter-free predictions of the NG model and of the Luscher term formula for the first 
excited states. We substitute the values of a 2 a that appear in Table fl] in Eq. (13.4p and 
Eq. (13. 5p for w — 1, q — and Nr = Nl = 1 and compare with our data. We present the 
comparison for SU(3) and SU(6), and for two different lattice spacings in Fig. [31 

In contrast to the case of the ground state, both the predictions Eq. (13.51) and Eq. (13. 6p 
are unable to describe our data for all the values of I that we study, and its seems that they 
will become consistent with the data only when I ~ (6 — 7)/ \fo ~ 2.7 — 3.2 fm. Despite this, 
the NG prediction is strikingly within ~ 5% of our data for I > 2/ \fa ~ 0.9 fm, and becomes 
consistent when I > 3.5 jyfa ~ 1.6 fm. To further check how this agreement depends on iV 
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FIG. 3: The energy of the 1st excited state, divided by al (See Tables HHTI1 for the values of 
the string tensions). Left panel: Comparing SU(3) and SU(6) for a similar lattice spacing of 
a? a ~ 0.03. Right panel: Comparing different lattice spacings for SU(3). The black line is the 
NG prediction, while the dotted(dashed) lines are the NG prediction expanded to leading and 
next-to-leading order. 

and whether it evolves with the lattice spacing a, we have repeated this analysis for gauge 
groups with 2 < iV < 8 and a large set of lattice spacings with 0.2 fm > a > 0.05 fm, and 
string lengths with I > 3/y/a, and this will be presented in a forthcoming publication 25]. 
What we find is that the dependence on both a and N is not significant. 



VI. COMPARISON OF THE NG SPECTRUM WITH THE LOWEST ~ 30 STATES 



In this section we present results from an extensive calculation of the lowest ~ 30 states 
in the spectrum. We restrict the discussion to states with w — 1, q — 0, 1, 2 and Nl,r < 3. 
We have obtained results for the third excited states with q = and N L = N R = 3 as 
well, but we postpone their presentation to j^. The results were obtained for SU(3) 
with (3 = 21.00,40.00 and for SU(6) with (3 = 90.00. Here we have focused on strings 
with 1.4/ ' \fo < / < 5.5/ y/a. All energies that we present here were obtained from single 
exponential fits. 

We present the results in Figs. HE The lines are the predictions of NG model. The string 
tensions used for these predictions were extracted from the ground state energies with the 
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1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 2 2.5 3 3.5 4 4.5 



FIG. 4: The energies of the lowest 7 states in units of \fa as a function of Zi/<7- The three lines are 
the NG predictions for the ground state(red), I s * excited state(blue) and 2 nd excited state(black). 
Left panel: We present the results for the case of SU(3) for two different values of (5 (two different 
lattice spacings). Right panel: Results for SU(Q) with (3 = 90.000. In both panels we denote the 
degeneracy of the states in the legends. 

use of the fitting ansatz Eq. ( 13. 8ft . and we present the fitting parameters in Table ILT1 



Gauge group 


(3 


a 2 a 


c 3 


Confidence level 


SU(3) 


21.00 


0.030258(26) 


0.160(21) 


69% 


40.00 


0.007577(13) 


0.05(31) 


15% 


SU(6) 


90.00 


0.029559(36) 


0.04(21) 


88% 



TABLE II: The parameters a 2 a and C3 in the fit Eq. (|3.8|) . which are obtained for the S data sets. 
The errors on C3 are much larger then the ones that appear in Table U because here we did not 
study the very short strings. 

It is clearly seen that the NG predictions are very good approximations to the flux-tube 
spectrum and deviate from our data only at the level of ~ 2% once / > 4.2/ y/a ~ 1.9 fm. 
At this level some of the systematic errors may be significant. Note that the degeneracy 
pattern predicted by the NG model is seen from our data : the second energy level is fourfold 
degenerate at large-/. This degeneracy includes two positive parity states and two negative 
parity states, and these start splitting significantly once I £ ~ 1.35 fm. 

We note in passing that the 1/7 and 1/7 3 contributions to the energies, which were pre- 
dicted using effective field-theories are clearly insufficient to describe our data, and would 
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presumably do so only for much longer strings. In contrast, the implicit resummation of 
powers of the Luscher term which appears in the covariant string expression (see discussion 
at the end of Section IIIIBI) is strongly supported by our data. 

Next, we fitted the data for the excited states. In the case of the first excited energy 
level, where there is only one state per level, we used the fitting ansatz Eq. (13. 8ft . In the 
case of the second excited energy level, where for each parity there are two states, we used 
a modified form of Eq. (13.81) and fitted the difference between the energies squared of these 
states to the form 

(SEf^j^. (6.1) 

The results of these fits are presented in Table IIHI 

To compare our data to the Luscher- Weisz prediction in 4j we momentarily assume that 
p = 2. Expanding Eq. (I3.8P this assumption results in an 0(1 /I 3 ) contribution to the string 
energy which deviates from the Luscher- Weisz prediction by an amount proportional to 
In the case of the first excited state and I > 2/y/a, we substitute the values C 2 — 3 — 8 
from Table and find that this deviation is only at the level of 2% — 6%. The smallness 
of this deviation, and the largeness of the errors we find for p from our fits (see Table IIIII ) , 
demonstrate that to unambiguously determine p, and test the Luscher- Weisz prediction, 
requires statistical errors which are at least 2-3 times smaller than the ones our current data 
has, and a simultaneous control of any systematic errors that may be important at this level 
of few percents accuracy. 

Finally, we move to the nonzero longitudinal momentum q sector. As mentioned in 
Section IIIII in the NG prediction of Eq. ( 13.41) . the number of left and right movers are 
constrained to obey the level matching condition Eq. (13.11) . The comparison of this prediction 
with our data is given in Fig. 0, where we present \J E 2 /a — (2iiq/ ^fal) 2 as a function of 
/•v/a. 4 As clearly seen the data is indeed very well described by Eq. (13. 4p . 

VII. CONCLUSIONS 

We have calculated the energy spectrum of closed strings in SU(N) gauge theories in 
2 + 1 dimensions with the string length I in the range 0.45 fm < / < 3 fm. 



This way of presenting the data separates the data sets and eases the representation. 
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Level 


gauge group 


P 


ly/a > 


P 


C p 


X 2 /d.o.f 




SU(3) 


21.00 


2.1 


1.8(5) 


6(3) 


2.4/7 


i excrGcti 


SU(3) 


21.00 


1.4 


3.7(2) 


37(5) 


39.5/9 
















SU(3) 


40.00 


2.1 


1.7(7) 


3(2) 


2.0/2 




SU(6) 


90.00 


2.1 


1.6(4) 


8(3) 


3.8/5 




SU(3) 


21.00 


2.1 


3.0(3) 


53(15) 


4.2/7 


2 nd excited, P = + 


SU(3) 


40.00 


2.1 


0.5(5) 


4(2) 


3.2/4 




SU(6) 


90.00 


2.1 


0.3(3) 


3(1) 


3.4/5 


z excited, r — — 


SU(3) 


21.00 


2.1 


2.1(7) 


16(9) 


7.9/7 
















SU(6) 


90.00 


2.1 


2.5(5) 


36(18) 


5.5/5 



TABLE III: The parameter C p in the fit Eq. fl33|) and x 2 /d.o.f, which were obtained for the S 
data sets (see Section Hi C ip . For the second excited level we fit both the positive parity (P = +) 
and negative parity (P = — ) states. 

For the ground state we have studied 2 < iV < 8 for lattice spacings 0.06 < a < 0.11 fm, 
depending on N, and saw that the Nambu-Goto (NG) free bosonic string model describes 
our data very well even at the relatively short distances of I ~ 0.6 — 0.7 fm (this was already 
noted in [5(). In particular, we find that the central charge c is in general consistent with 
1. For example, in SU(2) we see c = 0.9991(55) when we extract it already at / ~ 0.68 fm, 
while in SU(5) we find c = 1.004(23) at I ~ 0.82 fm. This provides unambiguous evidence 
that the closed flux-tube is in the bosonic string universality class. 

To study the excited string spectrum we constructed a basis of ~ 80 — 200 operators 
in each channel, and from a variational calculation we extracted the lowest ~ 30 states. 
These include states with positive and negative parity, as well as with nonzero momentum 
along the string direction. In general we find that the agreement with the NG prediction in 
Eq. (13.41) is very good, including the expected degeneracy pattern. 

This agreement is in striking contrast to what we find when we simply compare to the 
Luscher term as in Eq. (13. 5p or to the Luscher-Weisz prediction in Eq. (13.61) . While for the 
ground states these describe our data already at / ~ 1.35 fm, where they are indistinguishable 
from the full NG formula, for the excited states the situation is completely different. In 
particular, whereas the NG prediction works well for the first excited state already at I ~ 1.35 
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N R = 3, N L = 1, q = 2, w = 1 



2 2.5 3 3.5 4 4.5 



FIG. 5: y E 2 /a — (2T\qj\foV) as a function of Zy^ for the ground state (g.s.) and excited states 
(e.s.) of the non-zero q sector. The four lines present the NG prediction Eq. (|3.4p . For Nr = 
3, Nl = 1 and q = 2, the expected degeneracy is three, which agrees with our excited state data, 
i.e. two(one) states with positive(negative) parity in red(blue). 

fm, the predictions of Eqs. (13.51) and (13.61) are still very far from the data. 

This means that our results show unambiguously that the confining flux-tube can be 
described by a covariant string theory with small/moderate corrections, down to very short 
distance scales, and possibly at all distance scales at large- N. 

Acknowledgements 

AA acknowledges the support of the EC 6 th Framework Programme Research and Train- 
ing Network MRTN-CT-2004-503369. BB was supported by PPARC and thanks the Isaac 
Newton Institute for Mathematical Sciences at Cambridge, in which this work was com- 
pleted. The computations were performed on machines funded primarily by Oxford and 



19 



EPSRC. 



[1] M. Liischer, K. Symanzik and P. Weisz, Nucl. Phys. B 173, 365 (1980), M. Luscher, Nucl. 

Phys. B 180, 317 (1981). 
[2] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, Phys. Rept. 323, 183 



(2000) |arXiv:hep-th/9905111|. 



[3] J. Polchinski, arXiv:hep-th/9210045 



[4] M. Luscher and P. Weisz, JHEP 0407, 014 (2004) [arXiv:hep- th/0406205 1, JHEP 0207, 049 

(2002) [arXiv:hep-lat /0207003| . 
[5] B. Bringoltz and M. Teper, PoS LAT2006, 041 (2006) |arXiv:hep-lat/0610035| , Phys. Lett. 

B 645, 383 (2007) [arXiv:hep-th/0611286| . 
[6] D. Karabali, C. j. Kim and V. P. Nair, Phys. Lett. B 434, 103 (1998) [arXiv:hep-th/9804132] , 

and references within. 

[7] J. Kuti, PoS LAT2005, 001 (2006) [PoS JHW2005, 009 (2006)] |arXiv:hep-lat/0511023] , 



N. D. Hari Dass and P. Majumdar, JHEP 0610, 020 (2006) |arXiv:hep-lat/0608024|. 
M. Caselle, M. Hasenbusch and M. Panero, arXiv:0707.0055 [hep-lat], P. Giudice, F. Gliozzi 



and S. Lottini, JHEP 0701, 084 (2007) |arXi v:hep-th/0612131| . B. Brandt and P. Majumdar, 
poster presented at the XXV Lattice Symposium, July 2007. 
[8] J. Polchinski and A. Strominger, Phys. Rev. Lett. 67, 1681 (1991). 



[9] J. M. Drummond, |arXiv:hep-th/0411017HarXiv:hep -th/0608109. N. D. Hari Dass and P. Mat- 



lock, arXiv:hep-th/0606265, arXiv:hep-th/0611215 



[10] M. Teper, [arXiv:hep-lat /9804008 1. Phys. Rev. D59:014512, 1999. 

[11] B. Lucini and M. Teper, Phys. Rev. D 66, 097502 (2002) [arXiv:hep-lat /0206027] . 



[12] B. Lucini, M. Teper and U. Wenger, JHEP 0406, 012 (2004) |arXiv:hep-lat/0404008| 



[13] H. Meyer and M. Teper, JHEP 0412, 031 (2004) [arXiv:hep- lat/0411039 1. 

[14] K. G. Wilson, Closing remarks at the Abingdon/Rutherford Lattice Meeting, March 1981, 
K. Ishikawa, M. Teper and G. Schierholz, Phys. Lett. B 110, 399 (1982). B. Berg, A. Billoire 
and C. Rebbi, Annals Phys. 142, 185 (1982) [Addendum-ibid. 146, 470 (1983)]. K. Ishikawa, 
M. Teper and G. Schierholz, Phys. Lett. B 110, 399 (1982). M. Luscher and U. Wolff, Nucl. 
Phys. B 339, 222 (1990). 



20 



B. Bringoltz and M. Teper, in preperation, and PoS LAT2007, 291 (2007) [arXiv: 0708. 3447] . 
B. Bringoltz and M. Teper, PoS LAT2006, 041 (2006) [arXiv:hep-lat /0610035| . 
H. B. Meyer, Nucl. Phys. B 758, 204 (2006) [arXiv:hep-lat /0607015| . 

J. Liddle and M. Teper, hep-lat/ 0509082 and in preparation. K. Holland, |hep-lat/0509041| J. 



Liddle, D.Phil Thesis, Oxford 2006. 

B. Bringoltz and M. Teper, Phys. Rev. D 73, 014517 (2006) [arXiv:hep-laV0 508021 1. 

Y. Nambu, Phys. Rev. D 10, 4262 (1974). Phys. Lett. B 80, 372 (1979). J. F. Arvis, Phys. 

Lett. B 127, 106 (1983). 

J. Polchinski, Cambridge, UK: Univ. Pr. (1998) 402 p 
P. Olesen, Phys. Lett. B 160, 144 (1985). 



B. Lucini and M. Teper, Phys. Rev. D 66, 097502 (2002) |arXiv:hep-lat /0206027 1 . H. Meyer 

and M. Teper, JHEP 0412, 031 (2004) |arXiv:hep-lat/0411039| . 

N. D. Hari Dass and P. Majumdar, JHEP 0610, 020 (2006) [arXiv:hep-lat70 608024] ■ 

A. Athenodorou, B. Bringoltz and M. Teper, To be published. 



21 



